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ABSTRACT 

The present work analyses the dynamics of a suspension of heavy particles in shear flow. 
The magnitude of the particle inertia is given by the Stokes number St = mj/Gmia, which 
is the ratio of the viscous relaxation time of a particle t p = m/Girqa to the flow time 7~'. 
Here, m is the mass of the part icle, a is its size, p is the viscosity of the suspending fluid and 
7 is the shear rate. The ratio of the Stokes num ber to the Reynolds number, Re = Pfio 2 /p, 
is the density ratio Ppjpj- Of interest is to understand the separate roles of particle (St) 
and fluid (Re) inertia in the dynamics of suspensions. In this study we focus on heavy 
particles, p p /pj > 1, for which the Stokes number is finite, but the Reynolds number is 
sufficiently small for inertial forces in the fluid to be neglected; thus, the fluid motion is 
governed by the Stokes equations. On the other hand, the probability density governing 
the statistics of the suspended particles satisfies a Fokker-Planck equation that accounts 
for both configuration and momentum coordinates, the latter being essential for finite St. 
The solution of the Fokker-Planck equation is obtained to 0(St) via a Chapman-Enskog 
type-procedure, and the conditional velocity distribution so obtained is used to derive a 
configuration-space Smoluchowski equation with inertial corrections. The inertial effects are 
responsible for asymmetry in the relative trajectories of two spheres in shear flow, in cont rast 
to the well known symmetric structure in the absence of inertia. Finite St open trajectories 
in the plane of shear suffer a downward lateral displacement resulting from the inability 
of a particle of finite mass to follow the curvature of the zero-Stokes-number pathlines. In 
addition to the induced asymmetry, the 0(57) inertial perturbation dramatically alters the 
nature of the near-field trajectories. The stable closed orbits (for St = 0) in the plane 
of shear now spiral in, approaching particle-particle contact in the limit. All trajectories 
starting from an initial offset of 0(57 1/2 ) or less (which remain open for St — 0) also spiral 
in. The asymmetry of the trajectories leads to a non-Newtonian rheology and diffusive 
behavior. The latter because a given particle (moving along a finite St open trajectory) 
suffers a net displacement in the transverse direction after a single interaction. A sequence 
of such uncorrelated displacements leads to the particle executing a random walk. The 
inertial diffusivity tensor is anisotropic on account of differing strengths of interaction in the 
gradient and vorticity directions. Since the entire region (constituting an infinite area) of 
closed orbits in the plane of shear spirals onto contact for finite St , the latter represents a 
singular surface for the pair-distribution function. The exact form of the pair-distribution 
function at contact is still, however, indeterminate in the absence of non-hydrodynamic 
effects. It should also be noted that finite St non-rectilinear flows do not support a spatially 
uniform number density owing to the cross-streamline inertial migration of particles. 
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Problem & Motivation 
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For suspensions where inertial forces are important much less is 
known and our predictive ability is severely limited. 
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the effects of particle and fluid inertia. This can only be done in 
microgravity because of the gravitational settling that occurs when 
the particle and fluid densities do not match. 


Related work 
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Theory and simulation for finite St, zero Re suspensions in 
sedimentation and shear by D.L. Koch and coworkers. 
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Dynamic simulation for nonzero Re , St, and (j> in 
shear flow. 


g Analytical theory: Re = 0, 0(50 
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Integrate out the momentum variables using the method of multiple scales 
(Chapman-Enskog-type expansion) to derive the configuration-space 
Smoluchowski equation with inertia. 



Method of Multiple Scales 
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+ V x -(R- 1 -F°g) = StV x -(Vg ) 
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zero -Stokes 
trajectory 
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Inertial Dynamics Simulation 

Early results on the shear viscosity and first and second normal stress 
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Conclusions 
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Inertial Dynamics Simulations with an O (N In N) method have 
been initiated. 



